Fidelity of holonomic quantum computations in the case of random errors in the 

values of control parameters 
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We investigate the influence of random errors in external control parameters on the stability of 
holonomic quantum computation in the case of arbitrary loops and adiabatic connections. A simple 
expression is obtained for the case of small random uncorrelated errors. Due to universality of 
mathematical description our results are valid for any physical system which can be described in 
terms of holonomies. Theoretical results are confirmed by numerical simulations. 
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I. INTRODUCTION 

The idea of holonomic quantum computations (HQC) 
was mainly developed in the work [lj. In holonomic 
quantum computer quantum gates are implemented as 
holonomies (parallel transporters alongloops) in the 
space of external control parameters A M [H y [j Practi- 
cal implementations of such computational scheme 0, Q 
are based on the systems described by some degenerate 
hamiltonian which depends on the set of external control 
parameters A M 0, 1 



H{X) = U{\)H U\\) 



(1) 



where U{\) are some unitary operators. The holonomies 
are obtained in this case by adiabatically changing the 
values of the control parameters in such a way that the 
initial and the final set of parameters coincide (thus one 
passes along the loop in parametric space) . The states of 
the quantum register are supposed to be the eigenvectors 
of the hamiltonian corresponding to a single degenerate 
energy level PHHH. 

However practical implementations of HQC encounter 
many difficulties, which are caused either by decoherence 
processes or imperfect external control. While decoher- 
ence processes in quantum information processing is a 
well-studied subject, the influence of errors in the values 
of external control parameters on the stability of HQC 
is specific to this particular implementation of quantum 
computations and attracted less attention. Improved in- 
sensitivity of HQC to external perturbations is believed 
to be one of its main advantages, therefore it is impor- 
tant to study the stability of HQC in as general case as 
possible. Some works on the topic are |E IE 0> @> El- 
Geometric phase in open systems was studied, for in- 
stance, in J7|, where the geometric phase was shown to 



be insensitive to dephasing. The effect of fluctuations 
in the classical control parameters on the Berry phase 
of a spin 1/2 interacting with magnetic field was ana- 
lyzed in |8| ■ Numerical study of the stability of a partic- 
ular implementation of HQC with respect to parametric 
noise was recently presented in the paper Analytic 
expressions for the fidelity of HQC implementation of 
Hadamard gate were obtained in [9J, revealing improved 
robustness against perturbations up to the fourth order 
in the magnitude of errors. Non-abelian Stokes theorem 
was applied to study the stability of HQC in the case 
of arbitrary loops and adiabatic connections in the work 
but the expressions obtained in this work do not take 
into account randomness of errors. 

In this work we will investigate the stability of HQC 
with respect to random errors in the values of external 
control parameters in the case of arbitrary loops and ar- 
bitrary adiabatic connections. The paper is organized 
as follows: in this section we introduce some notations 
which will be convenient for further calculations, in the 
section[H]we follow the work [BJ and use the non-Abelian 
Stokes theorem to obtain the expression for the fidelity 
of HQC. The section ITTT1 contains the main result of this 
paper, namely the expression for the fidelity in the case 
of small random errors. Finally, in the section IIVI the- 
oretical results are compared with results of numerical 
simulations. 

We will assume that the Hilbert space of the states of 
the quantum register used in HQC is A^-dimensional. In 
order to describe the evolution of the state vector of the 
quantum register one needs to introduce the adiabatic 
connection An which has the following matrix structure 



Etna 



(K) (A) = 

\ / ran 



v (A)|-i^|^(A)) 



(2) 
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where \4>m{X)}, m = 1...N are linearly independent 
state vectors of the quantum register. 

An arbitrary quantum gate is represented as the holon- 
omy of the adiabatic connection which corresponds to 
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some loop 7 in the space of external control parameters 



V exp 



dA"i M (A) 



(3) 



where V is the path-ordering operator. By choosing an 
appropriate loop 7 one implements different unitary op- 
erators (quantum gates) acting on the state vector of the 
quantum register: \<p) = P 7 | (pi n ), where \<f>i„) is the initial 
state of the quantum register. It should be noted that 
due to the mathematical universality of holonomy all our 
results can be actually applied to any physical system 
which can be described in terms of holonomies (the most 
interesting example is probably Yang-Mills field) . 

One can also state that evolution of the state vector of 
the quantum register is described by the parallel trans- 
port equation, which is equivalent to J3J: 



d_ 

ds 



(4) 



where s is some parametrization of the path 7 and t m — 
4r- is the tangent vector. 

as G 

Let us introduce the density matrix /5(A) = 
\4>(X)) (4>(X)\ of the quantum register for further conve- 
nience. Parallel transport of the density matrix is then 
described by the following equation: 



ds 



(5) 



It is convenient now to use the generators T a , a = 
1 . . . TV 2 — 1 of SU (N) group, which are normalized as 

Tr (r a Tbj = 5 a b- We use latin indices to denote the ele- 
ments of su(N) Lie algebra and greek indices for tensors 
on the manifold of external control parameters. The gen- 
erators T a in fundamental representation build a basis in 
the space of traceless hermitan operators acting on N- 
dimensional Hilbert space, therefore we can decompose 
the adiabatic connection and the density matrix p as: 

p = N- X i + f aP a , A„ = A°J + T a A«, 
p a = Tr (pf a ) , Al = Tr (i M T a ) , A° = N^Tr i p (6) 

where we have taken into account that Tr p = 1 and 
assumed summation over double indices. Equation (J5J 
can be rewritten in terms of p a and A a ^ as: 



— n a — T^na Ab „c 
^ s P — T Ly bc^flP 



(7) 



where C% c are the structure constants of su(N) Lie alge- 



bra 



Tb, T c 



After omitting the indices of 
p a and introducing the adiabatic connection in the ad- 
joint representation of SU(N) group a c = C a bc A^, 



Aadj f 



-A a J*> we can also rewrite the equation Q as: 



ds 



p = T »A a *p 



(8) 



We will also use the notation ,0(7), or ^(7) to denote 
the decomposition coefficients (|HJ) or the density matrix 
obtained as a result of parallel transport along the loop 
7. According to the equation JHJ ^(7) can be represented 
as: 



Pil) = "P ex P 



(/ 



dX^Afi 



(9) 



where pi n is the initial density matrix of the quantum 
register: p m = \(j> in ){^ in \ = AT -1 / + f a pf n , p a m = 

(<t>in\Ta,\<l>in) ■ 



II. FIDELITY OF HOLONOMIC QUANTUM 
COMPUTER 



The values of external control parameters can always 
contain unpredictable errors, therefore it is necessary to 
know how these errors can change the results of quan- 
tum computations. The influence of errors in the val- 
ues of the external control parameters on the stability of 
quantum computations was studied, for instance, in the 
works 0, . In order to measure deviations of the results 
of quantum computations from the desired result we will 
use the quantity named the fidelity, which is convention- 
ally used to characterize the stability of quantum systems 
with respect to external perturbations [3, |a> LLlL LL2| • Ex- 
ponential decay of the fidelity means instability of quan- 
tum system • Following the works [Til IT^ we define 
the fidelity as: 



/ = Tr 03(7')P(7)) 



(10) 



where ^(7') is the density matrix of the quantum register 
obtained as a result of parallel transport along the loop 
with errors 7' and p(pf) is the density matrix which cor- 
responds to the desired operation. The physical meaning 
of the fidelity is the probability of obtaining the correct 
result of quantum computations. We will assume here 
that the initial points of the loops 7' and 7 coincide, as 
it is usually true in practical implementations 0, H • It 
is also convenient to express the fidelity in terms of the 
coefficients p a in the decomposition (JSJl: 



f = N- 1 +p a {i)p a { 1 ) 



(11) 



Using the expressions © and (|11|) and taking into ac- 
count the properties of path-ordered exponents we obtain 
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for the fidelity: 



f = N- 1 +plVe X p 



( 



N- 1 



\ am I a „a 



<i:r"A"r | ,;//; 



(12) 



\S-y 



where £7 is the loop obtained by travelling forward in 
the loop 7 and backward in the loop 7'. 

For the purposes of further analysis we will apply the 
non-Abelian Stokes theorem [T3I IT3 to the expression 
(|12[) . as was proposed in the work y5j- In contrast to the 
work |{| we use adjoint representation of su(N) algebra, 
thus we must first define several new objects. We define 
the curvature tensor F^ v of the adiabatic connection as: 



0A«(A) dA«(X) 



dX v 



C a bc A b (X)A:(X) (13) 



The shifted curvature tensor G" is defined as: 



/ A 



G° (A,A ,x) =Vexp 



d\»Af \ a b F b (X) (14) 



where Ao is the so-called "reference point" and \ is some 
path which connects the points Ao and A. The shifted 
curvature tensor is the curvature tensor parallel trans- 
ported to the reference point. We also introduce the 
shifted curvature tensor in the adjoint representation 



° bc U HV 



Now it is possible to apply the non-Abelian Stokes the- 
orem to the path-ordered exponent in l|12l) and express 
the integral over the loop £7 as the i nteg ral over the sur- 
face Ss-y spanned on the loop £7 



ic mtcgr 

Mm 



f = N- 1 +plVs Sy exp 



\s 6 



a a 
b riri 



(15) 



where do^ v is the element of area on S*5 7 and Vs s -, is the 
path-ordering operator for the surface Ss 1 . The initial 
point of the paths 7' and 7 is chosen as the reference 
point Ao in the definition of the shifted curvature tensor 
(|14fl . The surface Ss 7 can be chosen to be the surface 
of the minimal area spanned on the loop #7. The paths 
X are assumed to lie on the surface Ss-y- The choice 
of the paths x completely determines the path-ordering 
operator for the surface Vs Sl 0, • 

The distance between the paths 7' and 7 should be 
small, as otherwise HQC will lead to unpredictable re- 
sults and will be therefore useless. Let us denote the 



difference between the coordinates of the corresponding 
points on the loops 7' and 7 as <5A M . Such definition im- 
plies that <5A is the function of the points on the path 7. 
To the first order in SX^ we can write the element of area 
as da^ v ~ dX^SX^ . We retain only the first-order terms 
in 5X 11 in the argument of the path-ordered exponent in 
(|15p. However we do not expand the whole exponent to 
the terms of order 8X V , as the total effect of numerous 
small errors can be quite significant for large loops 7, 
which can be inferred from the results of our numerical 
simulations (see further) . After these assumptions we can 
write, using the properties of the path-ordering operator 

v Sh Mm 



'Pssy exp 



Pexp / dX»8X v {X)G<${X) 



(16) 



Such approximation is valid if the shifted curvature ten- 
sor does not change significantly at distances of order <5A 
and if SX is small compared to the loop size. Further- 
more, to the first order in dX^ the path x m (II 5|> can be 
assumed to coincide with the part of the path 7 which 
connects the initial point with the point where the shifted 
curvature tensor is calculated. 

Thus we obtain the following approximate expression 
for the fidelity: 



+ Pin 



V exp 



dX l 5X v Gfj 



bPv 



(17) 



III. THE FIDELITY IN THE CASE OF 
RANDOM UNCORRELATED ERRORS 

Errors in the values of external control parameters are 
unpredictable and essentially random, which makes the 
general expression (|17f) for the fidelity quite impractical. 
Indeed, it is necessary to know the errors SX^ exactly 
in order to calculate the fidelity in the form (|17fl . One 
would like instead to predict the fidelity basing on the 
statistical properties of errors. This was done numerically 
in some particular cases 6]. Here we will obtain a general 
expression in the limit of small and uncorrelated errors. 

As the fidelity is the probability of obtaining the cor- 
rect result of quantum computations, the correct expres- 
sion for the fidelity which takes into account the random 
character of errors is obtained by averaging the expres- 
sion (|17l) over all errors according to the multiplication 
law for conditional probability: 



- 1 + pI, 



■Pexp 



dX^SX^G 



fiV 



(18) 
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where by the line over the expression we denote averag- 
ing with respect to errors in the values of external control 
parameters. In this work we will make the following as- 
sumptions concerning the statistical properties of errors: 

1. Expectation values of errors are equal to zero: 
6X^ = 0. 

2. Errors are statistically independent in the parts of 
the loop 7 separated by sufficiently large distance, 
so that one can introduce some characteristic cor- 
relation length A c which has the meaning of the 
distance between uncorrelated errors and which de- 
termines how frequently the errors occur along the 
loop 7. 



3. Errors are small: \\G* || |£A^| A c < 1, where the 
norm||G-„||= / £ (G^) 2 . 



It should be noted that ||G°„|| = \\F£ V \\, and the small- 
ness condition can be also written as: 



I|f;jm^|a c «i (19) 

The expectation value of the path-ordered exponent in 
(|18f) can be calculated using the van-Kampen cumulant 
expansion p^l[lr|: 



Pexp I J dX^8X»G ad J 



= 7>cxp (f>/fc! f ... f dAf 



dK h ((G^* 1 (A 1 )...G°* fc (A fc )«JA^ ...5A"*» ] (20) 

7 7 



where ((. . .)) denote the so-called cumulants For 
instance the second-order cumulant is ((AB)) = AB — 
A B. The summand with k = 1 is equal to zero because 
the expectation value of errors is equal to zero. If the 



smallness condition (|19|) holds, the contribution of high- 
order terms is small compared to the second-order terms, 
and they can be omitted in the van-Kampen expansion 
(|2Up. which leads to the following expression: 



Pexp / d\»5\»G a J? «7>exp 1/2 / / dA^A^A«(Ai)^(A 2 )G^(Ai)G^(A 2 ) (21) 



As the correlation function (5A a (Ai)<5A^(A2) is not small 
only for close values of Ai and A2, that is only for close 
enough points on the loop 7, it is possible to integrate 
over Xi in (|21|l . The integrand is large enough only in 
the vicinity of the size A c of the point A2, thus the path- 
ordered exponent in l|21|l can be approximately calculated 
as: 



such integration can be regarded as the implicit definition 
of the correlation length A c : 



dA^dA^A^Ai^A^) 



7 7 



A c / d\»T v (\)6\ a (\)8\P(\) 



(23) 



Pexp / dX»5X»G ad J 



Pexp 1/2 A c / d\»T"5\ a 6\PG a ^Gl d * (22) 



where t v is the tangent vector of the loop 7. Note that 



After applying the Van-Kampen expansion the final 
approximate expression for the fidelity reads: 

f*N- 1 + 



-plVexp ( 1/2 X C J dX^6X»SX0G^G a ^ ) fr&(24) 

7 
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The expression (|24() shows that the stability of quan- 
tum computation essentially depends on: 

1. The duration and closeness of errors (the factor A c 
in (El) 

2. The shifted curvature tensor on the space of exter- 
nal parameters G"* 

3. The expectation value of the square of the magni- 
tude of errors 

4. The length of the contour 7 in the space of external 
control parameters 

The fidelity should in general decay with the length of 
the loop 7. Fidelity decay rate is proportional to the 
expectation value of the square of errors and the square 
of the shifted curvature tensor. 



IV. RESULTS OF NUMERICAL SIMULATIONS 

Numerical simulations were performed in order to ver- 
ify the theoretical results of the previous section. We 
used two-dimensional Hilbert space of the quantum reg- 
ister and worked with generators of SU(2) group. Nor- 
malized generators of SU (2) group are T a = a a / \/2, a = 
1, 2, 3, where a a are the Pauli matrices, and the structure 
constants of su(2) Lie algebra are C% b — e a b c /y2. Para- 
metric space was also assumed to be two-dimensional 
with parameters A 1 = x, A 2 = y. We have chosen a very 
simple form for the adiabatic connection, which however 
yields nontrivial coordinate dependence for the shifted 
curvature tensor: 

A x =& 1 , Ay = & 2 , Al = V2, A 2 y = V2 

Fly = -F^ x = 2V2 (25) 

with all other components being equal to zero. As both 
Al di and A" dj do not commute with i*!? * , the shifted 

y ^y 1 

curvature tensor is not constant. Such adiabatic connec- 
tion is, for example, approximately equal to the adiabatic 
connection for the optical HQC based on Kerr medium 
(described in the section IV of the work 17]) at small 
values of displacement parameters (x and y) and zero 
squeezing parameter (denoted as r\ in the work |17| h 
We used square loop 7 of the size L x L. The compo- 
nents of the initial density matrix of the quantum register 
were p u = 1, p 12 = P21 = P22 = 0. 

Errors were simulated in the following way: a set of 
N err points on the loop 7 was chosen, then each of these 
points was shifted independently from the others by a 
gaussian-distributed random vector. To obtain the loop 
7' these shifted points were connected by a cubic spline. 
Distribution of errors was assumed to be isotropic, i.e. 
5\ a S\P ~ 8 a p. Parallel transport equation JI) was solved 
using the fourth-order Runge-Kuttha method. The value 
of the fidelity was calculated using the definition (|10|l . 



Finally the fidelity was averaged over the implementa- 
tions of errors. Theoretical value of the fidelity H24|) was 
also calculated numerically. It should be noted that we 
have used not the correlation length A c but the number 
of error points N err on the loop as the parameter in our 
simulations. Correlation length can be estimated as the 
length of the loop divided by the number of error points 
N err : A c w 4LN~l. When calculating the fidelity as 
the function of the loop length, N err was increased pro- 
portionally to L, so that the correlation length remained 
constant. 

As can be inferred from the figures ^ 121 an d El the 

results of the simulations are in a good agreement with 
the expression (|24|l . Theoretical values of the fidelity ob- 
tained using the expression (|24|l are plotted with solid 
lines. It is worth noting that the expression (|24|) involves 
matrix exponents rather than just scalars, and thus the fi- 
delity decay in general differs from simple exponential. In 
order to emphasize this difference exponential fits of the 
simulation results were also plotted on the figures [21 and 
|31 In figure[2]the fitting function is 0.5 + 0.5 exp (— aS\ 2 ), 
which is a naive estimation for l|24|) if one treats the expo- 
nent as a scalar. Thus the dependence of the fidelity on 
perturbation strength differs from gaussian. Similarly, 
dotted line on the figure [2] corresponds to the function 
0.5 + 0.5 exp (— (3/N err ). Here the exponential approxi- 
mation is also not very good. 

Figures © and J2J) show that there are small discrep- 
ancies between theoretical results and the results of the 
simulations at large correlation lengths A c and at large 
error magnitudes, which means that the value of l|19|) is 
not sufficiently small in this region, and the van-Kampen 
expansion with second-order terms can not be applied 
to calculate the expectation value of path-ordered expo- 
nent in H21fl ■ Let us estimate the smallness parameter for 
the extremal values of perturbation strength and corre- 
lation length used in our simulations. The norm of the 

curvature tensor is ||-F°„|| = (F^ v ) 2 — 4. When 

V a, n, v 

SX = 0.3, L = 10, N err = 400, A c k, 0.1, the smallness 
parameter is \\F£ V \\ SX A c w 0.12 (figures □ and 0). When 
SX = 0.15, L = 10, N err = 100, A c » 0.4, the small- 
ness parameter is ||F^J| SX X c pa 0.24 (figure OJ- Thus 
the results of the simulations indicate that the expression 
(|24|l is valid even if the fidelity is not close to unity, that 
is when the total effect of numerous small errors is not 
small. The results of our simulations are in agreement 
with those in |(|. 



V. CONCLUSIONS 

In this paper we investigated the influence of random 
uncorrelated errors on the stability of holonomic quan- 
tum computations in the most general case, without ex- 
plicitly specifying the adiabatic connection or the loop in 
the space of external control parameters. A simple ex- 
pression for the fidelity was obtained, which is valid when 
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the condition \\F£ V \\ \SM\ X c < 1 holds, where ||F^|| is 
the norm of the curvature tensor of the adiabatic con- 
nection. This expression can be used to predict fidelity 
decay rates by specifying just a few parameters such as 
the average square of parameter deviations and the fre- 
quency of errors. It should be noted that though the fin al 
expression is valid for small magnitude of errors |£A^|, it 
is applicable even when the fidelity is not close to unity. 
Our theoretical results describe with a good precision the 
data obtained by numerical simulations, some discrepan- 
cies arising only when the condition ljT§)l does not hold. 



FIG. 1: The fidelity as the function of the loop size L for 
different magnitudes of errors at fixed correlation length (A c = 
0.1). 
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FIG. 2: The fidelity as the function of the magnitude of errors 
\J S A 2 at fixed loop size (L = 10) and correlation length (A c = 
0.1). 
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FIG. 3: The fidelity as the function of the number of er- 
ror points N er r at fixed loop size (L — 10) and perturbation 
strength (<5A = 0.15). 



